Spin and density excitations in the triangular-lattice t-J model with multiple-spin 

exchange interactions: 3 He on graphite 
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Using an exact diagonalization technique on small clusters, we study spin and density excitations 
of the triangular-lattice t-J model with multiple-spin exchange interactions, whereby we consider 
anomalous properties observed in the doped Mott region of the two-dimensional liquid 3 He adsorbed 
on a graphite surface. We find that the double-peak structure consistent with experiment appears 
in the calculated temperature dependence of the specific heat; the low-temperature sharp peak 
comes from the spin excitations reflecting the frustrated nature of the spin degrees of freedom 
and high-temperature broad peak comes from the density excitations extending over the entire 
band width. The clear separation in their energy scales is evident in the calculated spin and 
density excitation spectra. The calculated single-particle excitation spectra suggest the presence of 
fermionic quasiparticles dressed by the spin excitations, with an enhanced effective mass consistent 
with experiment. 

PACS numbers: 67.30.-n, 67.80.-s, 67.30.hr, 67.80.dm 



I. INTRODUCTION 

3 He atoms adsorbed on a graphite surface is known to 
be an ideal two-dimensional correlated spin- 1/2 fermion 
system. A solidified commensurate phase of 3 He atoms 
is stabilized at a 4/7 density of the underlying layer of 
4 He atoms due to the substrate potential corrugation 
and thus a triangular lattice of 3 He atoms is formed, 
which is a realization of a gapless quantum spin liq- 
uid (QSL)ji Theoretically, this 4/7 phase of spin-1/2 
3 He atoms has been studied by using the triangular- 
lattice Heiscnbcrg model with the multiple-spin exchange 
interactions .SiMiSiSii Importance of the density fluctua- 
tions has recently been pointed out as well4 

A finite amount of vacancies of 3 He atoms can be intro- 
duced into this 4/7 phase in a stable manner, where the 
vacancies can hop from site to site of the triangular lattice 
via quantum-mechanical tunneling motions even at ab- 
solute zero temperature. The presence of such vacancies, 
called the zero-point vacancies (ZPVs), was predicted a 
few decays agoi^^ Quite recently, the experimental evi- 
dence for the ZPVs has been reported in the monolayer 
of 3 He adsorbed on a surface of graphite preplated by 
a solid monolayer of 4 Hei n i 12 i.e., heat capacity mea- 
surements of the system show an anomalous coexistence 
of a magnetic round-peak near 1 mK and a broad peak 
at several tens mK that are associated with the ZPVs 
doped into the commensurate Mott-localized solid. The 
ZPVs are maintained up to the doping of almost 20% of 
the lattice sites, which we call the doped Mott region of 
monolayer 3 He. 

Theoretically, Fuseya and Ogatai^ have proposed the 
triangular-lattice t-J model with four-spin ring-exchange 
interactions as an effective model for the doped Mott re- 
gion of the system and obtained its ground-state phase 
diagram. The low-energy excitations of the model were 
also discussed. They have thereby argued that there is 



a new-type anomalous quantum-liquid phase, character- 
istic of the "spin-charge separation" , which may be rele- 
vant with the anomalous features observed in the doped 
Mott region of the monolayer 3 He adsorbed on a graphite 
surface. 

Motivated by such developments in the field, we study 
in this paper the triangular-lattice t-J model with the 
multiple-spin exchange interactions further. In particu- 
lar, we directly calculate the spin and density excitation 
spectra and single-particle spectra as well as the tem- 
perature dependence of the specific heat and uniform 
magnetic susceptibility by using an exact-diagonalization 
technique on small clusters. We thereby consider the 
anomalous properties observed in the doped Mott region 
of the two-dimensional liquid 3 He adsorbed on a graphite 
surface. 

We will thus demonstrate that the double-peak struc- 
ture actually appears in the temperature dependence of 
the specific heat, which is quantitatively consistent with 
experiment; the low-temperature sharp peak comes from 
the spin excitations and high-temperature broad peak 
comes from the density excitations. The spectral weight 
for the calculated spin excitation spectra is concentrated 
on a very low-energy region that scales with the exchange 
interactions, while that of the density excitations extends 
over an entire band width that scales with the hopping 
parameter of the vacancy. The clear separation between 
spin and density excitations in their energy scales is 
thus found. The accumulation of the low-energy spectral 
weight of the spin excitations comes from the frustrated 
nature of the spin degrees of freedom of the system; i.e., 
the ferromagnetic two-spin interactions compete with the 
antifcrromagnetic four-spin interactions on the geomet- 
rically frustrated triangular lattice. The single-particle 
excitation spectra suggest that the vacancies behave like 
fermionic quasiparticles dressed by the spin excitations, 
with the enhanced effective mass consistent with exper- 
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iment. Preliminary results of our work have been pre- 
sented in Refpii. 

This paper is organized as follows. In Sec. II, we 
present our model and method of calculation. In Sec. Ill, 
we present our results of calculations for the specific heat, 
magnetic susceptibility, spin and density excitation spec- 
tra, and single-particle excitation spectra. We compare 
our results with experiment in Sec. IV. We summarize 
our work in Sec. V. 
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FIG. 1: (Color online) (a) Noninteracting tight-binding band 
structure of the triangular lattice with the nearest-neighbor 
hopping parameter t. Horizontal line represents the Fermi 
energy at half filling, (b) Brillouin zone and Fermi surface at 
half filling. 



II. MODEL AND METHOD 

The triangular-lattice t-J model with the multiple-spin 
exchange interactions is defined by the Hamiltonian 

« = -*E (*U* + H - c -) + J E ( s < • s . - ! r 1 ) 

(ij),<r {ij) 
+Kj^(P i + P^)+R (Ps + Ps 1 ) (1) 

(ijkl) (ijklrnn) 

where ci a = Ci a (l — 7i< -o-) is the projected annihila- 
tion operator of a fermion ( 3 He atom) at site i and spin 
a (=TjI)) allowing no doubly occupied sites, Si is the 
spin-1/2 operator, and rii (= mi + n^) is the num- 
ber operator. The summation in the t-J part of the 
model is taken over all the nearest-neighbor pairs (ij) on 
the triangular lattice. P4 and Pq are the four-spin and 
six-spin exchange operators defined as P4 = PuPikPij 
and P 6 = PinPimPilPikPij, respectively, where P lj = 



(1 + <Ji ■ 0j)/2 with the Pauli spin matrix The sum- 
mation is taken over all the possible combinations of four 
nearest-neighbor sites (ijkl) for P4 and over all the equi- 
lateral hexagons (ijklrnn) for Pq. 

In this paper, we study the dynamical properties of 
the model under the introduction of vacancies in the 4/7 
commensurate solid phase of 3 He, i.e., removal of parti- 
cles ( 3 He atoms) or addition of ZPVs. We thus define 
the filling n of particles as n = N / L where N is the total 
number of particles and L is the total number of lattice 
sites in the system; in particular, n = 1 is referred to as 
"half filling", which corresponds to the 4/7 solid phase. 
The noninteracting band structure and Fermi surface at 
half filling of the tight-binding model with the nearest- 
neighbor hopping parameter t are shown in Fig. 1; we 
find no nesting features in the Fermi surface and no sin- 
gularities in the density of states for n < 1. 

The nearest-neighbor hopping parameter t and two- 
spin and four-spin exchange interaction parameters J and 
K have been estimated as follows: 13 t ~ 50 — 100 mK, 
—J ~ 1 — 10 mK, and K/\ J\ ~ 0.2. We note that, in this 
parameter region, the two-spin exchange term favors the 
ferromagnetic spin polarization (J < 0) but the four-spin 
exchange term gives the antiferromagnetic spin correla- 
tions between neighboring spins, and thus we have the 
situation where the strong frustration in the spin degrees 
of freedom of the system appears. Geometrical frustra- 
tion also appears on the triangular lattice when the in- 
teraction between spins is antiferromagnetic. We have 
examined the effects of the six-spin exchange interac- 
tion term on the ground state and excitation spectra and 
found that the effects are very small, in particular when 
doped with vacancies. This is because the Pq exchange 
interaction is easily cut by the presence of vacancies. We 
will therefore present the results at R = in this paper, 
the model of which we refer to as the t-J-K model. 

Throughout the paper, we use t = 1 as the unit of 
energy unless otherwise stated and we set H = fee = 1. 

We use the Lanczos exact-diagonalization technique on 
small clusters to calculate the ground state and excita- 
tion spectra of the model. In particular, we calculate the 
dynamical spin and density correlation functions defined, 
respectively, as 

g(qjttf ) l^ |^t __l__ s S ^ ) (2) 

and 

AT(q )W ) = -l3<* |4— , 1 v Tlqltto), (3) 

TT ^uj + 17] — (H — h ) 

where $0 an d E are the ground-state wave function and 
energy, respectively. and n q are the Fourier trans- 
forms of the spin and particle-number operators defined, 
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FIG. 2: (Color online) (a) 20-site and (b) 12-site clusters used 
for calculations, (c) Brillouin zone and available momenta of 
the 20-site cluster in the periodic boundary condition. 



respectively, as 



Aq " vl 



(4) 
(5) 



where Tj is the position of the lattice site i. 

We also calculate the single-particle excitation spec- 
trum defined as 



^(q,^)^^ (q, -w) + A + (q,w) 



(6) 



with the particle removal spectrum 

,-4 1 

7T 



A-(q, w ) = —cj^lc^ __^c qCT |vl/ ) (7) 

H uj + irj- (H- E ) 



and particle addition spectrum 

a + (<l«) = -^o|c q ^ + ^_^_^ o) at CT |^), (8) 

where r\ — ► +0, which is replaced by a small positive 
number in the actual calculations to give an artificial 
broadening of the spectra. We use a cluster of 20 sites 
with periodic boundary condition for these calculations 
(see Fig. 2), where the independent available momenta 
in the Brillouin zone, qo, • ■ • , qio, are also shown. They 
are at q^ = mh\ + ni>2 with bi = 2-7r/5(l, 1/V3) and 
b 2 = (0, 7r/\/3), where (m,n) are (0,0) for q , (1,0) for 
qi, (3,-1) for q 2 , (2,-1) for q 3 , (1,1) for q 4 , (0,2) for 
q 5 , (0,1) for q 6 , (2,0) for q 7 , (1,-1) for q 8 , (2,-2) for 
qg, and (1, —2) for qio. In the following, we will in partic- 
ular examine the cluster with two vacancies, i.e., n = 0.9. 

To calculate the temperature T dependence of the spe- 
cific heat C(T), magnetization M(T) under uniform mag- 
netic field h, and uniform magnetic susceptibility x(T) = 



\im.h^o{9M I dK)T, the Hamiltonian for a smaller-size 
cluster of 12 sites (see Fig. 2) is fully diagonalized to 
calculate the partition function . 15 ! 16 We add the Zeeman 
term —h ^ Sf to the Hamiltonian Eq. (1) when we cal- 
culate the magnetic response of the system. In the fol- 
lowing, we will in particular examine this cluster with 
one vacancy (n = 0.92) and two vacancies (n — 0.83). 



III. RESULTS OF CALCULATION 
A. Ground-state phase diagram 

The ground-state phase diagram in the parameter 
space of the present model at R — has been ob- 
tained by Fuseya and Ogata^ which we have also re- 
produced successfully. Here, we briefly review their re- 
sults. The phases obtained are as follows: the region 
of phase separation (phase-I), the region of Fermi liq- 
uid with strong spin fluctuations (phase-II), the region 
of new- type anomalous quantum liquid (phase-Ill) , and 
the region of fcrromagnctism (phase- IV), where we follow 
their notations of the phases. They have put special em- 
phasis on the phase-Ill, which has been argued to be the 
region of "spin-charge separation" and may be relevant 
with the anomalous features of the doped Mott region of 
the 3 He monolayer. 

In the following, we in particular examine the region 
of this new-type anomalous quantum liquid (phase-Ill) 
using the parameter values J = —0.3 and K = 0.06, 
which we compare with the results of other regions when 
necessary, i.e., the region of ferromagnetism (phase-IV) 
using J = —0.3 and K = and the region of Fermi liquid 
(phase-II) using J = -0.3 and K = 0.15 - 0.2. 



B. Specific heat and entropy 

The calculated results for the temperature dependence 
of the specific heat C(T) at n = 0.83 and 0.92 are shown 
in Fig. 3. We find that there appears a double-peak 
structure in C(T) at K = 0.06 (corresponding to the 
new-type anomalous quantum-liquid phase); i.e., a sharp 
peak at low temperatures and a very broad peak extend- 
ing over high temperatures. The double-peak structure is 
not clearly seen at K = and 0.15 — 0.2. We should note 
here that the specific heat coefficient 7, where C(T) = 7T 
at low temperatures, cannot be deduced from the present 
calculations since C (T) decays exponentially at low tem- 
peratures due to the discreteness of the energies of finite- 
size systems. 

We will show in Sec. Ill D that the low-temperature 
sharp peak comes from the excitation of the spin degrees 
of freedom of the system and the broad high-temperature 
peak comes from the excitations of the density degrees 
of freedom of the system. In other words, the width 
of the sharp low-energy peak scales with the exchange 
interactions between spins (a combination of J and K) 
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FIG. 3: (Color online) Calculated temperature dependence of 
the specific heat (per site) of the t-J-K model at J = —0.3. 
We use the cluster of 12 sites with two vacancies (n — 0.83) 
in the upper panel and with one vacancy (n = 0.92) in the 
lower panel. 



and the width of the broad high-energy peak scales with 
the hopping parameter t of the vacancy. 

It is interesting to note that the double-peak struc- 
ture in the specific heat C(T) due to the separation in 
their energy scales between spin and density degrees of 
freedom has previously been discussed in the context of 
the low-energy excitations in the Hubbard ladder systems 
with charge ordering instability although the latter is for 
the insulating systems with a charge gap»il In Sec. IV, 
we will compare the obtained double-peak structure with 
experimen t 11 ' 12 . We also calculate the temperature de- 
pendence of the entropy S(T) (not shown here), which 
will be compared with experiment^, also in Sec. IV. 



C. Uniform magnetic susceptibility 

The calculated results for the temperature dependence 
of the uniform magnetic susceptibility x(T) are shown in 
Fig. 4. We find that the temperature variation is strongly 
dependent on the value of K: (i) When < K < 0.04, 




FIG. 4: (Color online) Calculated temperature dependence 
of the uniform magnetic susceptibility of the t-J-K model at 
J = —0.3. We use the cluster of 12 sites with two vacancies: 
n — 0.83. Straight line corresponds to the Curie law, x(T) = 
C/T, where C is the Curie constant. 



x(T) is strongly enhanced in comparison with the Curie 
susceptibility x(T) = C /T, resulting in the ferromag- 
netic spin polarization at low temperatures, (ii) When 
K = 0.06 at which the frustration in the spin degrees 
of freedom is the largest, x(T) is slightly suppressed in 
comparison with the Curie law. Here, the ground state 
is highly degenerate due to the frustration of the spin 
degrees of freedom; in our cluster, the degeneracy is 15 
fold, (iii) When K > 0.07, x{T) is rapidly suppressed 
with decreasing temperatures. Here, the ground state of 
the system is spin singlet without degeneracy. 

It should be noted that these results come basically 
from the finite-size effects of small clusters. However, we 
may infer the intrinsic nature of the infinite-size system 
and its K dependence from the low-energy behavior un- 
der the magnetic field. The results for x(T) thus obtained 
are compared with experiment in Sec. IV. 



D. Spin and density excitation spectra 

The calculated results for the spin and density exci- 
tation spectra at K = 0.06 are shown in Fig. 5. The 
excitation spectra for the corresponding noninteracting 
infinite-size system, iVo(q, ui) — 2So(q, w )> are a l so shown 
for comparison. We find the following. The spectral 
weight for the spin excitations is concentrated on a very 
low-energy region of around u> < 1 . This reflects the pres- 
ence of a large number of nearly degenerate low-energy 
states coming from the frustrated nature of the spin de- 
grees of freedom. The low-energy spectral weight is ex- 
tended over the entire Brillouin zone, rather than special 
momenta, reflecting the spatially localized nature of the 



5 



1 






K .... 


% 

j i . j —i i-i — ■ i 


Jf 

■> 




% 

*r 1,1, 






J 










" J 


t 


i 




J 






ft' 1 





2 4 6 8 10 



FIG. 5: (Color online) Calculated spin (solid line) and density 
(dashed line) excitation spectra of the t-J-K model at J = 
-0.3 and K = 0.06. t] = 0.04 is assumed. We use the 20- 
site cluster with two vacancies (n — 0.9). The spectra for 
the noninteracting infinite system iVo(q, u>) are also shown 
for comparison (thin solid line). 



spin fluctuations. The spectral weight for the density ex- 
citations, on the other hand, extends over a wide energy 
range of about < u> < 9 (entire band width), which is 
more or less resembles the spectrum of the noninteracting 
system. 

From these results, we may say that the spin and 
density excitations are clearly separated in their energy 
scales: i.e., the spin excitations concentrate on the low 
energy regions, the width of which scales with the ex- 
change interactions (a combination of J and K), and 
the charge excitations extend over the entire band width, 
which scales with the hopping parameter t of the vacancy. 
With increasing K, we find the upward shift of the low- 
energy spectral weight of the spin excitations; e.g., at 
K = 0.15 — 0.2, we find the peaks at a higher-energy re- 
gion of around < uj < 5, where the momentum depen- 
dence of the positions of the peaks becomes significant as 



well. Thus, the separation between the energy scales of 
the spin and density excitations becomes weaker. 

We should note that the energy range where the spec- 
tral weight of the spin excitations accumulates, i.e., < 
w < 0.1, corresponds well to the temperature range where 
the low-temperature sharp peak in the calculated specific 
heat C(T) appears. We should also note that the broad 
spectra extending over the entire band width correspond 
well to the very broad high-temperature peak in the cal- 
culated result for C(T). We may therefore conclude that 
the separation between the spin and density excitations 
in their energy scales is responsible for the double-peak 
structure of the temperature dependence of the specific 
heat. 
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FIG. 6: (Color online) Calculated single-particle excitation 
spectra of the t-J-K model for the ground state of the 20- 
site cluster with two vacancies (n = 0.9). Vertical dotted 
line represents the Fermi energy. Vertical bars represent the 
position of the noninteracting band dispersion with a reduced 
hopping parameter t e g = t/6. We assume J = —0.3, K = 
0.06, and r\ = 0.04 



E. Single-particle excitation spectra 

The calculated results for the single-particle excitation 
spectra A(q,w) are shown in Fig. 6. From the results, 
we can deduce the possible quasiparticle band structure 
and hence the Fermi-surface topology. We find the fol- 
lowing. There are broad and incoherent spectral features 
over a wide energy range corresponding to the total band 
width of the noninteracting dispersion, i.e., —5 < u> < 5, 
but there emerge the sharp quasiparticle-like peaks with 
a characteristic dispersion in the vicinity of the Fermi 
energy. This result is similar to the case of the square- 
lattice t-J model near half filling j 18 i 19 Let us assume this 
to be the consequence of the presence of fermionic quasi- 
particles. Then, we find the quasiparticle band structure 
to be fitted well by the noninteracting band dispersion 
with a reduced hopping parameter t c s or with an en- 
hanced effective mass m* of the quasiparticle (see Fig. 7). 
From the fitting, we find the value 

7TL 

i cff = — t ~ (l/6)t (9) 

m 

or m* /m ~ 6 at n — 0.9 and K = 0.06. The present re- 
sults also suggests that the Fermi-surface topology of the 
quasiparticles is equivalent to that of the noninteracting 
system since the quasiparticle band can be obtained only 
by assuming that the band width is reduced (or the band 
mass is enhanced). Thus, the Fermi surface is large (i.e., 
its area oc n) rather than small (i.e., its area oc (1 — n)). 
The doping dependence of m* should be interesting, in 
particular whether m* diverges or not at n — ► 1. How- 
ever, we cannot answer this question in our small-cluster 
study; the behavior of m* even in the square-lattice t-J 
model still remains to be a puzzle. 

We may also assume that the quasiparticle band width 
shown in Fig. 7 may scale well with the energy of the spin 
excitations, i.e., a combination of the exchange param- 
eters J and K, while the entire band width (~ 9t) of 
the broad spectral features shown in Fig. 6 scales with t, 
as in the case of the square-lattice t-J model near half 
filling.— ^ Thus, we again find the separation between 
the spin and density degrees of freedom in their energy 
scales. The quasiholes (or quasiparticles as their con- 
jugate) are thus the vacancies dressed by the spin ex- 
citations. We should note therefore that the spin and 
density degrees of freedom are not exactly separated in 
this sense, unlike in the Tomonaga-Luttinger liquid in 
the one-dimensional interacting fermion systems^ Only 
the energy scales are different. Further experimental and 
theoretical studies will be required to clarify the true low- 
energy physics of the system. 

We may also point out that the enhanced effective mass 
of the quasiparticle band structure may partly be respon- 
sible for the enhancement of the effective mass m*/m 
determined from the specific heat coefficient 7 although 
the latter cannot be obtained from our finite-size calcula- 
tions. In Sec. IV, we compare the effective mass obtained 
from the quasiparticle band dispersion with experiment. 
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FIG. 7: (Color online) Quasiparticle band structure of the 
t-J-K model. Darkness of the shadow is in proportion to 
the spectral weight A(q, ui) shown in Fig. 6. Solid circles 
and bars represent the peak positions of the single-particle 
spectra. The peak positions are fitted with the noninteracting 
band structure (solid line) with a reduced hopping parameter 
i c ff = t/6. Horizontal line is the Fermi energy. 



IV. COMPARISON WITH EXPERIMENT 

The calculated results for the heat capacity C(T) at 
J = -0.3 and K = 0.06 (see Sec. Ill B) at the fillings 
of n = 0.92 and 0.83 are compared with experiment in 
Fig. 8(a). We here assume the value of t determined so 
as to reproduce the higher-temperature peak observed in 
C(T) at n = 0.92, i.e., t = 43.8 mK, so that we have 
I J| = 13.1 mK and K = 2.63 mK with keeping the ratio 
t : I J I : K — 1 : 0.3 : 0.06. We also assume that the total 
area of the sample used in experiment (556 m 2 ) is uni- 
formly active and contributes to the heat capacity. We 
should note that the results for the temperature region 
T < 1 mK are not reliable because of the finite-size ef- 
fects where the discreteness of the energies in the system 
gives the exponential decay of the heat capacity at low 
temperatures. 

We then find the fair agreement with experiment; in 
particular, the double-peak structure in C(T), i.e., the 
lower-temperature peak that comes from the spin excita- 
tions and higher-energy peak that comes from the density 
excitations of the system, are reasonably well reproduced. 
More precisely, we find that our calculated results repro- 
duce the experimental tendency that, near half filling, the 
lower-temperature peak is high but with increasing the 
vacancy concentration, the higher-temperature peak be- 
comes larger and simultaneously the peaks shift to higher 
temperatures. 

Note that the specific heat coefficient 7 (or the effective 
mass) cannot be estimated from the present calculations 
of C(T) due to finite-size effects. However, we find that 
the value of the enhanced effective mass m*/m ~ 6 esti- 
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FIG. 8: (Color online) Temperature dependence of (a) the 
heat capacity C(T) in units of mJ/K and (b) entropy S(T) 
divided by the entropy of N free spins Nk~B In 2. Comparisons 
between theory (solid and dashed lines) and experiment (sym- 
bols) are shown. Experimental data are taken from Rets . 11 ' 12 . 



The calculated results for the entropy S(T) are also 
compared with experimentally determined^ entropy in 
Fig. 8(b). We again find the fair agreement in their gen- 
eral tendencies. More precisely, we find that the cal- 
culated curves of C{T) cross the line of the entropy of 
N free spins 7Vfceln2 at ~10 mK, at which the lower- 
temperature peak in C(T) terminates. This result also 
supports that the lower-temperature peak in C (T) comes 
from the excitations of the spin degrees of freedom of 
the system. The higher-energy peak should therefore 
come from the density degrees of freedom or motions 
of vacancies in the system. Note that the experimen- 
tally determined 12 entropy is significantly smaller than 
2V&b1h2 even in the vicinity of half filling, n = 0.997, 
and even at temperatures of 10 — 20 mK where the lower- 
temperature peak in C(T) terminates. The missing en- 
tropy may reside in the region of much lower tempera- 
tures that the present experiment does not approach. 21 




r(mK) 

FIG. 9: (Color online) Temperature dependence of the mag- 
netization M under the uniform magnetic field h. Compar- 
ison between theory (solid lines) and experiment (symbols) 
is shown. The straight line corresponds to the Curie law. 
Experimental data are taken from Ref.— . 



mated from the calculated quasiparticle band structure 
(see Sec. Ill E) is consistent with the experimental value 
~ 7.5 — 10 estimated from the observed temperature de- 
pendence of C(T) at n = 0.89^ Here, we should note 
that the definition of m* in the experimental specific heat 
coefficient is two-fold: one is the value deduced from the 
lower-temperature peak and the other is the value de- 
duced from the higher-temperature peak. The two values 
are, however, not very different at least for n < 0.9, so 
that we can make comparison with our theoretical value. 
Then, it seems reasonable to assume that the renormal- 
ization of the band structure due to the spin excitations 
is mainly responsible for the observed enhancement of 
the effective massji 2 - Thus, stated differently, the specific 
heat coefficient 7 should be determined predominantly 
by the spin excitations of the system. 



The calculated results for the magnetization under the 
uniform magnetic field are compared with experiment in 
Fig. 9. We should note that the experimentally applied 
magnetic field, which is h = 0.006 in our calculations, is 
very small in comparison with the energy scales of the 
3 He system, so that the behavior of the magnetization 
under the uniform magnetic field is the same as that of 
the uniform magnetic susceptibility denned at h — > 0. 

We find that, although the finite-size effect is strong 
at low temperatures, the calculated magnetization at 
K = 0.06 is consistent with experiment in the sense that 
the value is somewhat smaller than the value expected 
from the Curie law M (T) = x(T)ft = Ch/T. For more 
quantitative comparison, however, one would need the 
techniques appropriate for treating infinite-size systems, 
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so that experimentally observed plateau-like behavior 22 
in the temperature dependence of the magnetization can 
be explained. 

V. SUMMARY 

We have used an exact-diagonalization technique on 
small clusters to study the low-energy physics of the 
triangular-lattice t-J model with the multiple-spin ex- 
change interactions, whereby we have considered the 
anomalous properties observed in the doped Mott region 
of the two-dimensional liquid 3 Hc adsorbed on a graphite 
surface. We have calculated the temperature dependence 
of the specific heat, entropy, and uniform magnetic sus- 
ceptibility, as well as the spin and density excitation spec- 
tra and single-particle spectra for the model, and have 
considered their implications. 

We have shown the following: 

(1) The double-peak structure appears in the temper- 
ature dependence of the specific heat. The result is 
quantitatively consistent with experiment. The low- 
temperature sharp peak comes from the spin excitations 
and high-temperature broad peak comes from the density 
excitations. 

(2) The spectral weight for the spin excitations is concen- 
trated on a very low-energy region, the width of which 
scales with the exchange interactions, while that of the 
density excitations extends over an entire band width, 
which scales with the hopping parameter of the vacancy. 
The clear separation between spin and density excitations 
in their energy scales is thus found. 

(3) The accumulation of the spectral weight of the spin 
excitations comes from the frustrated nature of the spin 
degrees of freedom of the system; i.e., the ferromagnetic 
two-spin interactions J compete with the antiferromag- 



netic four-spin interactions K on the geometrically frus- 
trated triangular lattice. 

(4) The single-particle excitation spectra suggest that the 
vacancies behave like the fermionic quasiparticles dressed 
by the spin excitations, of which the effective band mass 
is estimated to be m*/m ~ 6 at n — 0.9, in consistent 
with the effective mass measured from the specific heat 
coefficient. 

(5) The temperature dependence of the spin susceptibil- 
ity shows a suppressed Curie-like behavior, reflecting the 
situation where the ground state is highly degenerate due 
to the frustrated nature of the spin degrees of freedom. 

We hope that the present study will shed more light 
on the physics of the two-dimensional 3 He systems and 
stimulate further experimental and theoretical studies of 
the systems in greater details. We have focused on the 
doped Mott region of the monolayer 3 He in this paper. 
However, it has recently been reported^ 3 , that the bilayer 
3 He systems also contain rich physics concerning heavy 
fcrmions with quantum criticality, which we want to leave 
for future study. 
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